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Intensity interferometry with more than two detectors? 



m 

o 

(N 

Oh 
< 



6 

CO 

c3 



m 
m 

o 
m 

> 

X 



Vinay Malvimat, 1 Olaf Wucknitz 2 ' 3 and Prasenjit Saha 4 

1 School of Physics, University of Hyderabad, Gachibowli, Hyderabad-500046, India 

2 Max- Planck- Institut fur Radioastronomie, Auf dera Hugel 69, 53121 Bonn, Germany 

3 Argelander- Institut fur Astronomie, Auf dera Hugel 71, 53121 Bonn, Germany 

4 Institute for Theoretical Physics, University of Zurich, Winterthurerstr 190, 8057 Zurich, Switzerland 



12 April 2013 



ABSTRACT 

The original intensity interferometers were instruments built in the 1950s and 60s 
by Hanbury Brown and collaborators, achieving milli-arcsec resolutions in visible 
light without optical-quality mirrors. They exploited a then-novel physical effect, now 
known as HBT correlation after the experiments of Hanbury Brown and Twiss, and 
nowadays considered fundamental in quantum optics. Now a new generation of inten- 
sity interferometers is being designed, raising the possibility of measuring intensity 
correlations with three or more detectors. 

Quantum optics predicts some interesting features in higher-order HBT. One is 
that HBT correlation increases combinatorially with the number of detectors. Signal 
to noise considerations suggest, that many-detector HBT correlations would be mea- 
surable for bright masers, but very difficult for thermal sources. But the more modest 
three-detector HBT correlation seems measurable for bright stars, and would provide 
image information (namely the bispectrum) not present in standard HBT. 
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1 INTRODUCTION 

If two detectors are counting photons from an ordinary incoherent light source, there is a tendency for photons to arrive at 
both detectors together. This is known as photon bunching or HBT correlation, after the pioneering work of Hanbury Brown 
and Twiss, and is a consequence of the bosonic quantum statistics of photons. In a situation where coherent light would 
produce interference fringes, the HBT correlation with incoherent light varies according to those would-be but absent fringe 
patterns. This allows a type of interferometry with incoherent light and without optical-quality mirrors, known as intensity 
interferometry. 

HBT correlation appears naturally in classical wave optics, and was actually first used to build a radio intensity in- 



terferometer, which resolved an extra-galactic radio source for the first time (Hanbury Brown et al. 19521. But when the 



same ideas were applied to visible light, and the effect was measured, first in the lab (Hanbury Brown & Twiss 1956) and 
then with starlight (H anbury Brown fe Twiss|1958 l, it became controversial, because it implied that different photons could 
interfere, contrary to conventional wisdom at the time. The controversies were eventually resolved with the development of 



a quantum-statistical theory for incoherent light by Sudarshan ( 1963 1 and Glauber ( 1963 1 and the emergence of quantum 



optics. In quantum optics, ordinary light (or 'chaotic light') behaves like a random mixture of lasers, and the semi-classical 
picture of interfering waves, the squared amplitudes of which determine the emission rate of photo-electrons, turns out to be 
validrjThe general phenomenon of photon bunching or HBT correlation then spread to different areas of physics in different 
guises. For example, it is well known in the context of nuclear collisions ( Baym|1998 l, and it may even be relevant to animal 
vision ( |Sim et al.|2012 |. 

Meanwhile, Hanbury Brown and collaborators developed the Narrabri Stellar Intensity Interferometer (NSII) which 
measured stellar diameters down to milli-arcsecs (|Hanbury Browri]|l968). But the photon detectors then available were only 



1 The semi-classical description is not valid in general. A nice counter-example is provided by electrons, which show HBT anii-correlation 
(Kiesel et al. 12002 I from fermion statistics. 
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blue-sensitive, limiting the instrument to hot stars. So the NSII ran out of stars to observe in a few years, after which it 
was dismantled and almost forgotten. Decades later, with a new generation of detectors, intensity interferometry has become 



somewhat topical in astronomy again. A number of proposals and experiments have appeared in recent years (Ofir & Ribak 



2006a|b Borra|2008 Jain fc Ralston|2008 Foellmi 2009 Horch fc Camarata|2012 1. Particularly ambitious are plans to adapt 
Cherenkov telescopes into a giant reincarnated NSII with resolution down to < 0.03 mas (Nunez et al. 2012 Dravins et al. 



2012) 



Future intensity interferometers are likely to have many detectors, not just two. Clearly, detectors can work in simultaneous 

The basic 



pairs. The prospect of measuring JV-fold photon coincidences has also been suggested (e.g., Ofir fc Ribak| | 

theory is well established in the quantum-optics literature and three-detector correlations have been measured in laboratory 



experiments (Zhou et al. 20101. But are the results interesting for an astronomical instrument? This paper examines the 
question. 



2 STANDARD HBT 



A nice introduction to HBT from a modern quantum-optics perspective appears in Glauber ( 2006 1. We quote some key results 
here. 

The central concept is the field correlation functions. Let x\ and x 2 denote the spacetime locations of two detectors. The 
first-order correlation"] is defined as the average 

G(x 1 ,x 2 ) = (E ( - ) (x 1 )E (+ \x 2 ))^ (1) 

Classically, _B (±) denote the positive and negative frequency parts of the electric field at a detector. In quantum statistics, the 
fields become operators. We recognize G(xi, x 2 ) as the unnormalized correlation function or 'visibility' in radio interferometry, 
or the spatial Fourier transform of the source in the sky. Note that G(x 2 ,xi) is the complex conjugate (Hermitian conjugate 
in quantum statistics) of G(xi,x 2 ). Thus G(a;i,a;i), which is the count rate at xi, is automatically real and non-negative. For 
a laser source, G(xi,x 2 ) would be independent of time, but for chaotic sources, the correlation falls away over a coherence 
time 

Ar^l/Au. (2) 

Now consider a form of second-order correlation 

G <2) (xi,x 2 ,x 2 ,x 1 ) = {E ( -\x 1 )E < -'>(x 2 )E (+) (x 2 )E <+ \x 1 )) (3) 

which is the probability of coincident detection within At. Chaotic sources have the property (cf. Glauber's equation 37) 

G l - 2) {x 1 ,x 2 ,x 2 ,x±) — G(x 1 ,x 1 )G(x 2 ,x 2 ) + |G(xi, :r 2 )| 2 . (4) 

The first term corresponds to the random coincidences that would also be expected for classical particles, but the second term 
describes the non-classical HBT correlations. 

To estimate the signal to noise, suppose we have two detectors with equal count rates r (photons per time unit), 
parametrized as 

rAT = G(x 1 ,x 1 ) = G(x 2 ,x 2 ), (5) 

so that G(x, x) provides the number of photons per coherence time. Let us now count photons over some time At. The time 
resolution At of the detector, typically At 3> At, is often called the reciprocal electrical bandwidth, because in the early 
experiments it was set by amplifier properties. 

The product of photon counts (corresponding to random coincidences) will be close to (rAi) 2 , as given by the first term 
on the right of equation Q. The HBT signal given by the second term will make a small but non-zero contribution. Over 
a coherence time it will be |G(a;i, £2)| 2 - Over the much longer interval At there are, so to speak, At /At coherent slices, 
making the HBT signal r 2 Ar At. Using Poisson noise for the photon numbers and thus interpreting the square root of the 
total coincidence rate as the noise, we get a signal to noise ratio 

SNR(Ai)~rAr (6) 

for full correlation. This applies to a single counting time At. Over many counting times, the SNR adds in quadrature so that 
we get 

SNR(T)~tAt^ (7) 

when integrating over a duration T. When increasing the coherence time by using narrow-band filters, the photon rate r will 



2 In a different terminology our 'first-order correlation' would be called 'two-point correlation', and our 'second-order correlation' is a 
'four-point correlation'. 



© 0000 RAS, MNRAS 000, 000-000 



Intensity interferometry with more than two detectors? 3 

decrease with Av « 1/Ar so that the 'spectral intensity' rAr stays constant. This has the remarkable effect that adding 
narrow-band filters does not reduce the achieved SNR. This makes it possible to increase the SNR further by using spectral 
detectors that can distinguish between different photon energies and effectively record many narrow channels simultaneously. 
Or one can decrease the optical bandwidth so much that photon count rates are sufficiently low not to be affected by detector 
dead times. 

In the NSII, the counting time At was ~ 10 nanosec. Current off-the-shelf instruments can achieve At ~ 50 picosec (not 
to mention better quantum efficiency and broader wavelength response). Clearly the time resolution is very good, but it is 
still orders of magnitude shorter than the coherence time. Both these inequalities are important. 

The fact that At ^> At is what makes intensity interferometry interesting in the first place. In standard interferometry, 
optical paths have to be kept under control to better than A, which is extremely demanding mechanically. For intensity inter- 
ferometry, path differences do not matter as long as they are well below cAt, corresponding to metres in the first experiments. 
This not only relaxes the mechanical tolerances, it makes the signal immune to atmospheric fluctuations. Compared to other 
techniques, intensity interferometry is possible with simpler technology. Put in another way, if we lack the precision to measure 
the first-order correlation G(xi,X2) directly, we can still get information on it indirectly through the second-order correlation. 

The shortness of At, on the other hand, is what makes intensity interferometry usable with very low coincidence rates. The 
NSII could operate at r At ~ I0~ 5 , because it had 10 8 counting times per second and hence could build up SNR ~ 10 4 r At 
in a second. Current technology could deliver at least another order of magnitude better. 

The above description, in terms of counting photons, is well suited to optical astronomy. In radio astronomy, on the 
other hand, a description in terms of waves and intensities is standard. For incoherent sources, the _E (±) fields are considered 
as (complex) Gaussian random variables with autocorrelation functions determined by the characteristics of the receiving 
system. Each measurement naturally consists of a finite sum of random amplitudes, so that the fields and intensities vary with 
time even for sources of constant luminosity. This 'wave-noise' or 'self-noise' defines the fluctuations whose correlations are 
measured as the HBT effect. The wave-noise also adds to the photon shot-noise that is important in the optical domain of low 
photon rates, and this contribution actually dominates for high rAr. In radio observations the noise is generally dominated 
by a third contribution, the 'thermal' or 'receiver noise', which corresponds to the wave- noise of the signals produced by the 
receiving system itself. 

Here a conceptual difficulty arises: in radio astronomy, the electric field is considered as a classical field which can be 
measured, whereas in quantum optics the electric field is a non-Hermitian operator and hence not itself an observable. How 
to reconcile the radio-astronomy and quantum-optics pictures? One possible resolution is given in |Burke| ( |1969[ ). Here we 
suggest another, which goes as follows. Let there be a source field 5 (±) , and let us superpose it on a known local field L (±) . 
The resulting field 

£ <±, (x 1 ) = S (±) +L (±) (8) 

then gets its intensity measured: 

G(xi,xi) = (5 ( - ] S (+) ) + (L ( - } L <+) ) + <S ( -) L (+) ) + (£<"> S (+) ) _ (9) 

Here the first two terms on the right are just the intensities \S\ 2 and \L\ 2 of the two fields. The last two terms give a beating 
oscillation in the photon count rates, and it is from these oscillation that the phase S is inferred. Thus, even if the source 
field cannot be measured directly, it can be inferred indirectly. The SNR will then be the last two terms in Q divided by the 
noise. If \L\ ^> \S\, the noise in \L\ will dominate. If that noise is oc \L\, the SNR will be oc \S\. In other words, the SNR will 
be proportional to the square root of the source field intensity, just as expected from the photon-counting picture. 



3 MORE THAN TWO DETECTORS 

For chaotic sources, the higher order correlation functions are all given in terms of the first-order correlations. The result is 
presented in equation (10.27) of |Glauber| ( |1963[ ), which, with a slight modification of notation, reads 

N 

G (,l) (xi, ...,xn,xn,---,xi) — ^2 II G(x k ,Vx k ) . (10) 

V fc=l 

Here Vxk denotes the fe th element of a permutation of {xi, . . . , xn}- The sum is over all permutations. In each term, the first 
argument always runs as xi,x%,... whereas the second argument runs as a permutation of that ordering. As we remarked 
earlier, in quantum optics, the G are correlations between the field operators, but the semi-classical approach of treating 
the fields as classical and then interpreting intensities as photon probabilities is valid for light (cf. Sudarshan 1963). Indeed, 



equation ( 10 1 also appears in the classical theory of random Gaussian variables, where it is known as Isserlis' theorem. 
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For three detectors, the formula ( 10 \ gives 

G < - 3) (x 1 ,x2,x 3 ,x 3 ,x2,x 1 ) = G(xi,xi)G(x2,X2)G(x3,x 3 ) + G(x 1 ,x 3 )G(x2,x 2 )G(x s ,xi) + 

G(xi,X 2 )G(x2,X 3 )G(x3,Xi) + G(xi, xi) G(x 2 , x 3 ) G(x 3 , x 2 ) + (11) 

G(x 1 ,x 3 )G{x2,x 1 )G(x 3 ,X2) + G(xi,x 2 )G(x2,x 1 )G(x 3 ,x 3 ) . 
In each term here, the first argument always runs as x^,X2,x 3 whereas the second argument runs as a permutation of that 
ordering. One can rewrite this expression in another way, which is easier to interpret, by introducing the normalized correlations 

_ G(xi,x 2 )G(x 2 ,xi) 

G(xi,xi)G(x 2 ,x 2 ) ' . . 

_ G(x 1 ,x 2 )G(x 2 ,x 3 )G(x 3 , zi) 
G(x 1 ,x 1 )G(x2,x 2 )G(x 3 ,x 3 ) 
The three-point coincidence rate is then 

1 + 012 + 913 + ff23 + 23? 3123 (13) 

times the chance coincidence rate. Here the gij terms are the (normalized) power spectrumr] while the last term is the bispec- 
trum, whose phase is well known in radio astronomy as the closure phase. Thus, whereas two-detector intensity interferometry 
only measures amplitudes but no phases, multi-detector combinations are actually sensitive to certain combinations of phases 
and provide qualitatively new information. 

If the detectors are close together, all the terms are equal (py = g±2 3 = 1), resulting in a six- fold enhancement over the 
chance coincidence rate. This has been measured in lab experiments (Zhou et al.|2010l. At this point, one may get the idea 



from the JV! terms in ( |10| , that splitting up a single collecting area into N parts will increase the count rate. But in fact that 
will not happen. Splitting reduces the count rate in each detector by a factor of N, hence the JV-point coincidence rate would 
be N\ (r At/N) n oc y/~Ne~ N (r At) , using Stirling's approximation for large N. 

Let us now estimate the SNR. We have to be careful here, because the combinatorial formulas define the coincidences 
over At. We are interested in coincidences over At, and we have seen before in the case of two detectors that random and 
HBT coincidences scale differently. To take care of this, we have to understand that the intensities that are being correlated 
are integrated over a duration At consisting of many intervals of At, and correlations exist only for fields within the same 
short interval. The iV-point HBT signal will therefore be 

S i,..jvX(rArf(At/Ar). (14) 



We should emphasize that ( 14 1 is not the number of coincidences, but the number that remains after subtracting off the 
chance coincidence rate and all the lower-order HBT effects. Meanwhile, the chance coincidence rate is 

(rAt) N . (15) 

Hence 

SNR(iV, At) ~gi... N x(r Ar) N/2 {Ar/dtf 2 - 1 (16) 

For N = 2 and g\2 — 1 we recover the simple expression (IM for standard HBT. 

As we mentioned above, the iV-detector coincidence rate includes all the lower-order HBT signals as well. So, in principle, 
iV-detector coincidences could be used to extract the two-point HBT signal. Is it advantageous to do so, compared to standard 
HBT? To answer this question, let us consider the number of iV-detector coincidences over At 

r N At N +r N At N - 1 ArJ29 1 k + ... (17) 

j<k 

The random coincidences are given by the first term, and two-point HBT gives the next term. In the same way as before, we 

may derive 

SNR(iV -> 2, At) ~ T x (rAt) N/2 ^ , (18) 



where T denotes the sum in equation (171. Taking many counting times in quadrature and assuming full correlation {gjh = 1), 
we have, analogous to equation ¥j\, 

SNR(iV - 2, T) ~ iV ^ (r At)-/ 2 ^ ^ ■ (19) 

If we keep the total collecting area constant but split it into TV detectors, so that the individual counting rate goes with r/N 
we have 

SNR(7v^2,T) /.AA^- 1 

SNR(2,T) [ ' \ N J { ' 

3 In terms of normalized visibilities 7^, we have 312 = I712I 2 and 3123 = 712723731- 
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In the photon-counting regime, rAt/N < 1, and hence splitting detectors does not help. The situation is different in the 
case of rAt/N > 1, in which we could not count individual photons anymore but can still measure intensities. The SNR 
would then apparently increase with N. Formally, the optimal number of detectors for rAt = 1/10/100 is N = 3/7/41 with 
SNR benefits compared to two detectors of 1.15/14.6/(1.42 x 10 9 ). These numbers look extremely promising. But there is a 



paradox: equation (20 1 has SNR increasing with At (that is, with coarser time resolution), when rAt is large enough. Clearly, 
this estimate based on photon shot-noise cannot represent the whole truth. In fact, for very high photon rates, a new source 
of noise enters: it is the well-known wave noise from radio-astronomy, and as mentioned in the previous section, it can be 
considered as HBT correlation at a single space point over different times. We leave a full calculation of the case of high 
photon rates for a future paper, but the basic conclusion is that SNR over one counting time At cannot exceed unity, no 
matter how bright the source or how ridiculously large the light collectors. 



4 FLUXES AND COUNT RATES 

From the previous section, we see that the number of photons detectable in a coherence time, or r At, is central to the SNR. 
In a passband around v, a source with flux density F v gives 

rAr~— A, (21) 

hv 

where A is the collecting area. As mentioned before, the effect of Av cancels between r and At. Rewriting the expression, 

mixing wavelength and frequency, in the form 

' A "" u »"(^)(i)(ff) <22 > 

makes it easy to compare sources. 

Bright stars can be F„ ~ 10 4 Jy, but A ~ 10 -6 m. Hence r At <C 1 per m 2 . More detailed estimates are shown in Figurelll 
On the other hand, bright masersjcan have F v ~ 100 Jy at A ~ 1 cm, so a large dish easily receives r At > 1. Now, intensity 
interferometry is not needed to resolve masers, because at radio wavelengths it is much easier to control the delays sufficiently 
accurately for standard interferometry. But high-order HBT using astrophysical masers would be a novel quantum-optics 
experiment. 

Blackbody sources lead to a rather elegant estimate of r At. Consider a blackbody source at temperature T and angular 
area Q on the sky. We detect photons from it using a collector of area A, in a passband AA. For convenience, we will work in 
terms of a logarithmic passband 

A(lnA) = ^- (23) 

A 

We also write 

"m- <*> 

Recall the number density of a photon gas 

4^A(lnA) 

This is for one polarization state; since we are making rough estimates here, we will disregard the second polarization state. 
To get the photon flux, we multiply this by c£l /(4n) or XuQ/(4tt). The photon arrival rate in an area A is thus 

We can rewrite A A in terms of the coherence time q2p . The number of photons received in a coherence time is then 

We can also write r At in another way. Consider the baseline needed to resolve the source. The square of the baseline is 
the area v4 a i ry ~ A 2 /SI of an aperture whose Airy disc corresponds to the size of the source. This is a natural limit for the 
collecting area of detectors to just resolve the source and, e.g., measure its size, because larger detectors would require longer 
baselines and reduce the correlations g so much that the source would be 'resolved out'. With this definition, we have 

^^t^^T' (28) 



4 For natural masers to show HBT correlation, it is essential that they are not single-mode systems like artificial masers, but rather like 
chaotic superpositions of laboratory masers. 
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Figure 1. Photon flux (above) and photon flux per coherence time (below) for three blackbody sources, corresponding roughly to Sirius, 
Arcturus and Betelgeuse. The lower panel leads to simple estimates of signal to noise. 
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We see that for an HBT baseline adapted to the source size (that is, A/A^y fixed), the SNR only depends on AT. Thus, 
by going far into the Rayleigh- Jeans tail, rAr can in principle be made arbitrarily large (with the consequence of increasing 



wave noise), but in the brightest part of the spectrum rAr <C 1. Equation (28 1 also tells us that intercontinental baselines 
are not an option for HBT, because then A/A a i Ty would inevitably be very very small, driving down the SNR with it. 

A simple physical interpretation for the expression ( 28 1 is obtained by noting that it is A/A a i ly times the phase-space 



density of photons in blackbody radiation. As result of Liouville's theorem, this phase-space density is conserved when the 
radiation leaves the sources to travel towards the observer. When increasing the distance, the total flux gets diluted, but at 
the same time the apparent size of the source shrinks so that the momentum-space density increases, which compensates for 
the former effect. A a iry is a reciprocal measure for the apparent size of the source, and A/A^y by definition takes care of 
both effects. With increasing distance, the light bucket has to grow to pick up more photons, but can also form a smaller field 
of view and be susceptible to a smaller part of momentum-space. 

Figure [l] shows the count rate and rAr for three blackbody sources. These have T = 9940 K and angular diameter 
9 = .007", w hich approximates Sirius, T = 43 00 K, 9 = .02", similar to Arcturus, and T = 3500 K, and 9 = .04", similar 
to Betelgeuse. |Hanbury Brown fc Twiss (19581 measuring Sirius had a collecting area of about 2 m 2 , and quantum efficiency 



15% at .4 micron. From Figure n] we would predict SNR ~ 3 x 10 5 per counting cycle, less absorption and other losses. The 
reported value is SNR = 8.5 in 345 minutes, using 5-45 MHz, which corresponds to SNR ~ 1 x 10 -J per counting cycle. 
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5 DISCUSSION 



We may expand the title of this paper to two questions. First, is the quantum-optical effect of three-point or higher-order 
HBT measurable for any astronomical source? If so, would it tell us new things about the source? 



Measuring three-point HBT for bright stars appears feasible. Using equation ( 16 1 for N — 2 and N = 3 scaled to an 
integration time T, 

SNR(T, N = 2) ~ 5 i2 x (r At) \ 

— ( 30 ) 

SNR(T,iV = 3) ~ pi23 x (rAr) 3/2 ^ • 

For visible light with a narrowband filter, the coherence time At ~ 10 _12 s. Off-the-shelf photon counters can reach a time 
resolution of At ~ 10~ 10 s. From Figure 111 a square metre of collecting area gives rAr ~ 10 -4 from a bright star. These 



numbers suggest SNR ~ 1 in an hour. For quicker results one would want to increase the collecting area — from (161 we see 
that T oc A~ N for a given SNR. As always in HBT, optical-path tolerances need only be <C cAt. 

Three-point HBT would provide the three-point closure phase, or bispectrum, of the source on the sky. Is that worth 
having? Since two-point HBT gives only the power spectrum, and no phase information, having the bispectrum is likely to 
be an important advantage for image reconstruction of bright stars. 

Going to four or more detectors, for bright stars or any other thermal sources, would be difficult. The geometric estimate 



(28 1 of the photon count per coherence time indicates that r At <C 1 for any thermal source, except far in the Rayleigh- Jeans 
tail. But r At appears at progressively higher powers in the SNR. For non-thermal sources, the situation may be very different. 
In particular r At > 1 appears achievable for bright masers. Masers can be imaged using standard radio-telescopes, so HBT 
may not provide any new information on them. Nonetheless, it would be an interesting physics experiment to look for the 
combinatorial enhancement of the HBT effect for large N. The total number of photons would not increase, of course, they 
would just get more and more bunched. 
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